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We study how periodicity manifestations recently found in the turbulent globally coupled maps 

depend on the global feature of the couplings. We examine three non-locally coupled map models. 

In the first two, the all-to-all interaction is maintained but the coupling decreases with distance in a 

power and an exponential law. In the third, the interaction is uniform but cut off sharply. We find 

" " " that, in all three and in dimension D = 1,2,3, periodicity manifests universally from turbulence 

^— «, . when the same suppression of the local mean field fluctuation is achieved by the non-local averaging. 
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O 1. Introduction 

:z;. 

^— ^ ' The globally coupled map lattice (GCML) is one of the basic models of complex systems with interacting chaotic 

^— I , elements. The simplest GCML consists of N identical logistic maps with uniform couplings via their mean field h. It 

has only two parameters, the non-linearity a of element map and the coupling e of the averaging interaction. Yet it 

exhibits a rich variety of phases under the balance between the randomness and the coherence |I| and provides us with 

a simple testing ground of basic features of wide variety of coupled random elements such as the Josephson junction 

>«• I array, vortex dynamics, multimode lasers as well as biological networks. In particular, in its 'turbulent reginre' — 

(~| ' the region of high non-linearity and very weak coupling — the system at large N exhibits intriguing properties in 

'^ I evolution. 

ti ' Firstly, the maps are under unfailing weak coherence even if the parameters are so chosen that no visible syn- 

I chronization of maps occurs [g|. As a result the mean square deviation (MSD) of the time series of mean field ht 

^^ • is sizably enhanced from the value predicted by the law of large numbers {Sh'^ oc 1/N), while the ht distribution is 

^ ' Gaussian following the central limit theorem 0. This was called as a hidden coherence and triggered much scrutiny 

OO |§f§[- Furthermore, it has been recently found that the periodic windows of the element map systematically foliate 

and produce various periodicity manifestations (PM's) M^- Therefore, GCML in the turbulent regime is a system 

which sensitively mirrors the periodic windows with the background hidden coherence. The amazing fact that even at 

f*^ very weak coupling the system easily form periodic cluster attractors may have important consequences in a complex 

T-H ' system of coupled chaotic elements, especially in the activity of the brain. 

^D ' In this Letter we study to what extent the PM's depend on the all to all uniform coupling feature of GCML. We 

r^ will show that they occur universally in three non-locally coupled map models and their strength can be predicted by 

•^ I a simple estimation of the fluctuation of local mean fields around the overall mean field. 

^ ■ The PM's are organized by the maps by synchronization when the parameters a, e are in a balance that allows a re- 

'^ I duction of the high A^-dimensional dynamics to that of the element logistic map in a periodic window at non-linearity 

• 'H b. The balance defines foliation curves on the (a, e) plane and all GCML on the curves are universally governed by the 

KS window dynamics. The most prominent PM's are induced by the period three window 0. If GCML is on the foliation 

^ . curves of the p3 window, the maps organize themselves into almost equally populated three clusters, which oscillate 

■ ■ ' mutually in period three — p3c3 maximally symmetric cluster attractor (MSCA) M. With slightly higher e at the 

same a, that is, on the foliation curves from the intermittent region, the maps organize themselves in p3c2 cluster 

attractor. The p3 attractors are formed at any reduction factor r {r = b/a < 1). Generally, for a small reduction 

(r ^ 0.95), maps form p = c MSCA and p > c clusters respectively along the curves of period p window and in the 

nearby higher e. In MSCA, the MSD of the mean field fluctuation is minimized due to the high population symmetry 

and all observed MSCA's are linearly stable [0. Contrarily, in p > c states, the MSD turns out extremely high due 

to missing clusters to fulfill the orbits. For a large reduction (r < 0.95), the clusters are no longer formed but their 

renmants induce the same structure in the MSD — a valley and peak respectively along the curves of a window and 

in the nearby higher e. We use below the MSD curve as a representation of PM's. 

2. Non-locally coupled map models 

The GCML on the lattice A is 

xpit + 1) = (1 - s)fixpit)) + sht, P e A (1) 
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with the mean field ht = (1/-^) X^oeA /(^q(0) ^^'^ /(■'-) ~ ^ ^ '^^'^ ■ This is a two step process; the independent 
mapping followed by an interaction via the mean field ht with an overall coupling e. By adding (nl) over P we find a 
relation -^ X^pga ^'p ~ W X^PeA fi^p) — ^^^ mean field is kept invariant in the interaction. All the non-local models 
we consider below respect this invariance rule. 

2.1. A power law model: POWa 

As an extension of GCML let us consider a model 

x'p = (1 -e)f{xp) +ehp, 

hp=Y^ WpQf{XQ) 

QeA 

where each map couples to other maps via a local mean field hp. The Ap(P) is a set of maps at an equal distance p 
from a site P. For simple analytic estimates below, we approximate it by a set of points on the boundary of a {2p+l)^ 
square (cube) for D — 2(3). The number of maps in Ap, Up, is then 2, 8p, 24p^ + 2 respectively for D — 1,2, 3. We 
impose the periodic boundary condition and the maximum 'radius' Pmax of Ap is {N^/^ — l)/2. As a requisite the 
weights for a local mean field hp must add to one. This constrains the coefficients c^") and d^"^ as 



Pn 
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+ d(")5(") = 1, with 5(°) = Y^ uplp"-. (3) 



Under this constraint, one obtains with a simple arithmetic both the above invariance rule and a further important 
relation -^ X^PeA ^p ^ ^ — ^^^ average of the local mean fields is nothing but the mean field of the whole system. 

Let us make (||) into a model which interpolates the GCML and the nearest neighbor CML. In order to match with 
GCML at a = 0, the coefficient must be c'-^-' = d'-"^ = 1/N. In order to match with the nearest neighbor CML 

x'p = f{xp) + ^—l J2 /(■^q)-«i/(^p)| (4) 

at a -^ oo, the coefficients must be 0^°°-' — d^°°^ = l/{ni + 1). In both limits, c = d. Therefore, we set c^"-' = d*^"^ for 
all a as the simplest interpolation. Normalizing the couplings by (0), we obtain a one parameter extension of GCML, 
POWq, with hp given by 

V P=^ QeA,(p) / 

As shown in Fig. [^ PM's in POWa diminish with the increase of a. We choose eight marking points (a-h) and in 
the insets we display corresponding MSD curves. By the strength of PM's the a interval may be divided into three 
regions. (I) — One can observe all PM's that occur in GCML. From the GCML limit {a — 0) up to the point a, the 
full strength PM's are produced. The relevant dominant windows are marked on the MSD curve at a, which agrees 
precisely with the curve in GCML m. From a to d, the peak valley structures, except for that due to p3 clusters, 
gradually diminish. The peak due to p5 window starts diminishing at a and it becomes half-height at b. At c all the 
sub-dominant peak- valley structures vanish, and even the p5 peak vanishes at d. (II) — The region of p3 PM's only. 
It starts from d and the p3c2 peak disappears at e. At f, only a broad MSD peak is seen in the MSD curve. (Ill) — 
Essentially the region of the hidden coherence. Only broad MSD peak can be seen around the foliation zone of the 
p3 window. At the start of III and near the top of the peak, the correlator of maps decreases exponentially in time in 
a p3 motion. At g, the broad MSD enhancement becomes half-height and the correlator fails to sense the periodicity 
everywhere. At h, the MSD enhancement disappears. The transition points Ti, T2, T3 between the regions are d, f, 
h respectively. As we see in the bar chart, the variation of PM's occurs most quickly in Z? = 1 and it is prolonged in 
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FIG. 1. The bar chart for the PM's in POWq with a = 1.90 and the MSD curves at the eight marking points (a-h). I: full 
PM's, II: only the p = 3 PM's, III: only the hidden coherence. 



higher dimensions. We note that ax^ « 0.9, 1.9, 2.9, approximately in the ratios 1 : 2 : 3 for D = 1, 2, 3 respectively. 

2.2. A coupled map lattice with exponentially decaying couplings: EXPp^ 
Similarly, we obtain a model with hp given by 

1 



1 + S'(Po) 



f{xp) + 51 ^P^P« XI -^(^q) ' 



(6) 



p=i 



QeAp 



where Wp,po = exp(— (p — 1)/ po) is the exponentially decaying coupling and S^^"^ — X]p=r "'p'w^p,po- This reduces to 
GCML at po ~* oo and the nearest neighbor CML at po — * 0. The PM's diminish with decreasing po via the same 
process as above. 

2.3. A coupled map lattice with an interaction range k: CML^ 

Above two models maintain all-to-all coupling feature of GCML. Let us now consider a non-local CML with hp 
given by 



p=iQeAp(P) 



(7) 



Here K = {2k + 1)^ is the number of maps within a range k. We find that the PM's diminish with decreasing k, again 
via the same patterns of MSD curves. Furthermore, we find that remarkably the same PM's occur irrespective to the 
dimensions if the neighborhood encloses the same number of maps. For instance, the range k at Ti is 77—92, 5 — 6, 2 — 3 
in Z? = 1,2,3 respectively, but the number of maps K within k is 155 — 185, 121 — 169, 125 — 343 in Z? = 1,2,3. 
The large error in D = 3 comes from the large-step increase of K with k. Thus we determine the marking points in 
CMLk by a refined neighborhood; a set of lattice points Q around P with X]i=i(^Pi)^ < k^. In Fig. g we show the 
three regions by a bar chart in terms of K. We find the bars in 13 = 1,2,3 agree each other well. CML^ was used in 
an analysis of hidden coherence from the view of 'beat' of mean field H. The Fourier power spectrum of ht is shown 
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FIG. 2. Right: Bar chart of PM's in CML„ at a = 1.90. I: full PM's, II: only the p = 3 PM's, III: only the hidden coherence. 
The error bars compare the ambiguities in determining the transition points by two types of neighborhood. Center: The MSD 
curves at the pointed positions. Left: Power spectrum of the mean field at (a, e) = (1.99, 0.10) in _D = 1. 



also in Fig. 0. Interestingly, the Fourier peaks due to the beat become outstanding in accord with the onset of PM's. 



3. The universality of PM's in non-local models 



The difference between GCML and other non-local models is only in the interaction step. In GCML, the maps 
contract uniformly to /i by a factor 1 — £, while in others a map f{xp) is contracted to the local mean field hp 
which distributes around the overall system mean field h. Therefore, when the variance of ^p = hp — h is large, 
some distortion of map configuration must unavoidably be introduced in the interaction step. Contrarily, when the 
variance is small, such a distortion will be avoided and the non-local system may evolve just in the same way with 
GCML. Thus, it is natural to consider that the deviation from the global limit is controlled by the variance of £,p. 

The ^p is an weighted sum of maps of the form 



^p = E (^PQ - ^) /(^q)- 



q«ea 



(8) 



where Wpq is the couplings in (||), (g), (Q). If the spatial correlation between the maps are negligible, the variance 
of ^p niay be estimated at each time t as 



(C|,)a = {{hp - h)')A « T{{fp - h)')A 



T^Y.'^WpqY 



N' 



(9) 



where (• ■ •)a denotes the average over the lattice and J^q Wpq = 1 is used. The factor T represents the suppression of 
the variance of the ^p by taking the weighted mean of map values over the lattice A. At the global limit, Wpq -^ 1/N 









TABLE I. The leading A^ estimate of the f in POWc. 
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and ^ — > (strictly no variance). For intermediate couplings and large N the factor 1/N may be neglected and !F is 
solely determined by the couplings. Combined with the above consideration let us propose a working hypothesis that 
PM's occur universally in all non-local models when the factor T is the same and put it under scrutiny below. 

In CMLk, the suppression factor ^ is simply JF — 1/K — 1/N. This succinctly explains the observation that PM's 
occur with the same strength at common K in all D and uniformly diminish with decreasing (increasing) K {T). 

In POWq, the JT is given by 

The leading A^ estimates for T^' are tabulated in Table. |. We find in particular J^ « logiV/4A^ at a = 1/2, 1, 3/2 
for D = 1,2,3 respectively. This gives a prediction that the PM's would be universal among POW^^-^^2: POW^T^^ 

and P0Wfj3%. 

This is indeed the case; the full strength PM's, the same with those in GCML, are realized in all the three. The 
J-j^" in EXPpQ may be obtained by substituting S^" and lUp „^ to S^^ and 1/p^" respectively. 

In Fig. we compare POWq. (EXPp^) with CML^ in the right (left). The inset illustrates the case of the transition 
point Ti(d) in _D = 1 as an example. The curve is J- for POWq in dlQ). The measured a at Ti gives the vertical band 
taking account for the ambiguity in judging the MSD curve pattern. Hence, the crossing of the curve and the vertical 
band gives the estimate of IF in POW^ at its Ti in Z? = 1. On the other hand, T is universal over D in CML^. The T 
at Ti of CMLk gives the horizontal band, again counting for the ambiguity and averaged over D. Thus, the vertical 
axis is used for both J^'s, that in POWq and that in CML^. If both models share exactly the same J^ at Ti, the 
curve will pass through the crossing junction of the horizontal and vertical bands. In this example, the curve crosses 
the vertical band slightly above the junction and the estimated !F are (8.5 ± 1.5) x 10^'^ and (5.5 ± 0.5) x 10"^ in 
POWq and CML^.; respectively. Or, one can predict the a at Ti in POWq from K at Ti in CML^ using the J- curve 
of POWq. The prediction is 0.83 ± 0.02 to be compared with the measured 0.90 ± 0.03. In the overall comparison, 
only the junctions are shown by error-bars. We find that the hypothesis remarkably works with respect to all the 
marking points and in D = 1, 2, 3 for T ranging from 10"'' to 0(1). 

A few remarks are in order. 

(i) The curves in POWq agree approximately with each other after a scale transformation 1 : 1/2 : 1/3 in a up to 
Tu ss 2 X 10~^. This extends the rule obtained by the leading TV calculation. The same strength PM's occur up to 
the marking point e if a is in the ratio of the system dimensionality, just like the universal PM's at the same K in 
CML,. 

(ii) The horizontal bands exhibit three transition points observed in CML^ using refined neighborhood. In other 
models with the coarse neighborhood, the T3 is missing in Z? = 2 and both T2 and T3 are missing in Z? = 3. (See 
Fig. |l| for POWq). The curves of T explain the difference succinctly; they are constrained by the limiting values 1/3^ 
so they can pass through only the lowest two (one) bands in £) = 2(3). We have numerically checked that the missing 
transition points are retrieved in both POWq and EXPpg with the refined neighbors. 

(iii) Our estimate of IF is based on an approximation that the spatial correlation is negligible. This is a legitimate 
approximation. Firstly, we note that the even at the formation of cluster attractors such as p3c3 MSCA and p3c2 
states, the spatial distribution of maps does not show any visible clustering. To avoid a confusion we stress that the 
clustering of maps is in the map values and not in the spatial distribution. We have checked that over the whole 
turbulent regime of the three models no visible spatial clusters are formed. Furthermore we have numerically checked 
that the estimate (|^) works remarkably well in three models. For Z? = 1 — 3 and over e = — 0.1 the measured ratio 
(Cp)a/(('5/p)^)a plotted versus F both in double logarithmic scale lies on a straight line for F from 10~^ to 0(1), 
which fully covers the range of F in Fig. 0. There turns out some spread of data points only in the e region for the 
p3c2 cluster attractor. This is due to the variation of the ratio of map populations in two clusters when formed from 
different initial configurations. 

4. Conclusion 

In this note we have focused our attention to the recently found periodicity manifestations in the turbulent regime of 
GCML. We have conducted an extensive statistical analysis in three non-locally coupled map lattices over Z? = 1, 2, 3 
and examined to what extent they depend on the non-locality of the models. We have noted that the essential 
deviation of the non-local CML from the GCML stems in the variance of the local mean field around the overall 



w 



10 



1 - 



10 



10" 



10 




lO' lOP 10^ 10^ lO' 



CML 



Po 



GCML 



a 



6 8 
- CML 



FIG. 3. The comparison between POWa and CML„ (right) and EXPpg and CML^ (left) with respect to eight marked 
changes (a-h) of the PM's. For D = 1,2,3 respectively the J-d is shown by the solid, dashed, and dotted curves and 
A^ = 51^,51^,13^. The inset illustrates the method for the case of d (Ti) in D = 1. The a (po) axis is in the normal 
(logarithmic) scale. 



mean field. We have analytically estimated the suppression factor J-' of the variance under an approximation that 
the spatial correlation of maps in the turbulence regime is negligible and checked that this JF remarkably agrees with 
the numerical result. We have found a salient universally that, irrespective of the difference in construction and the 
dimension of the lattice, the periodicity manifestations occur at the same strength to a good approximation once JF 
is the same. 
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